Impact of scale-dependent bias and nonlinear evolution on the ISW 
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I summarize recent results from Smitii, Hernandez-Monteagudo & Seljak[ll], a study of tlie 
^\ • impact of nonlinear evolution of gravitational potentials in the LCDM model on the Integrated 

' Sachs- Wolfe (ISW) contribution to the cross-power spectrum of the CMB and a set of biased 

■ tracers of the mass. We use a large ensemble of A'^-body simulations to directly follow the 
potentials and compare the results to analytic perturbation theory (PT) methods. The PT 
predictions match our results to high precision for k < 0.2/iMpc~^. We analyze the CMB- 
density tracer cross-spectrum using simulations and renormalized bias PT, and find good 
agreement. The usual assumption is that nonlinear evolution enhances the growth of structure 

, and counteracts the linear ISW on small scales, leading to a change in sign of the CMB-LSS 

■ cross-spectrum at small scales. However, PT analysis suggests that this trend reverses at late 
times when the logarithmic growth rate f = d\n D /d\na < 1/2 or Q,rn{z) < 0.3. Numerical 

. results confirm these expectations and we find nonlinear enhancement of the ISW signal on 

' small scales at late-times. On computing the total contribution to the angular spectrum, we 

I find that nonlinearity and scale dependence of the bias are unable to infiuence the signal-to- 

• ■ noise of the current and future measurements. 
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The temperature fluctuations in the Cosmic Microwave Background (CMB) are directly sensitive 
to the presence of Dark Energy, through the change in energy that a CMB photon experiences as 
it propagates through an inhomogeneous Universe with time evolving gravitational potentials, 
', There are three main processes that give rise to changing potentials: the Integrated Sachs- Wolfe 

I effect jl, ISW], which deals with the linear evolution of potentials; the Rees-Sciama effect [1, 

RS], which is concerned with the nonlinear evolution of potentials; and the Birkinshaw-Gull 
effect 0], where time evolving potentials arise due to mass flows. We define the "nonlinear ISW 
Q"^ I effect" as the sum of all three contributions. 

Unfortunately, the nonlinear ISW effect can not be observed directly in the CMB auto-power 
spectrum, owing to the fact that it affects the low multipoles, where cosmic variance is large. It 
^ ' can however be directly observed by cross-correlating the CMB with tracers of the Large-Scale 

^ . Structure (hereafter LSS). [5i]. This analysis has recently been performed by a number of groups 

using the WMAP data and several LSS measurements (e.g. SDSS, NVSS, 2MASS), with claims 
of up to 4(7 level detections @, In a recent paper, [1] measured the cross-correlation between 
super-structures and super-voids with the CMB. On stacking the signal they found a ~ 4.5a" 
detection, in multiple WMAP bands, and the sign of which appeared consistent with late-time 
ISW. This result is very puzzling when one considers signal-to-noise (S/M) calculations within 
the LCDM model. These show that the maximum S/J\f ~ 7 for full sky surveys [5(]- How then 
is it possible to obtain such high S/Af, given the partial sky coverage of current surveys? One 
proposed solution is that nonlinear evolution of potentials and galaxy bias may increase the 
S/M. 

In this study we calculate the impact of nonlinear evolution of gravitational potentials on the 
angular cross-spectrum between a biased set of density tracers and the CMB. We ask whether 
such effects influence S/M. We pursue a two-pronged attack: our first line of inquiry is analytic, 
and we use perturbation theory (PT) to predict and this we do in ^ our second line is to 
use a large ensemble of A^'-body simulations to directly follow <!>, and this we do in f|3l 



2 The ISW effect in nonlinear PT 



The ISW effect may be written as jl] : 
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- / dm{h,x;t) , (1) 

where fi is a unit direction vector on the sphere, <I> is the dimensionless metric perturbation in 
the Newtonian gauge, which reduces to the usual gravitational potential on small scales, the 
'over dot' denotes a partial derivative with respect to the coordinate time t from the FLRW 
metric, x is the comoving radial geodesic distance X = / cdt/a{t), and so may equivalently 
parameterize time, to Us denote the time at which the photons are received and emitted 
(i.e. last scattering), respectively, c is the speed of light and a{t) is the scale factor. 

The rate of change of ^ can be calculated from Poisson's equation (V^<^ = AnGpSa'^): 
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Thus we require knowledge of the time evolution of 6 and its growth rate 5. 

The collapse of cosmic structures can be followed into the nonlinear regime using perturba- 
tion theory (PT). The solutions for S in Fourier space may be written 

<5(k,a) = f;[D(a)]"5„(k,ao) ; 5„(k) ^ | niklg^i^^5l^(k-qi-. . .-q„)F^^^ 

(3) 

where 5i(qj) represents an initial field at wavenumber qj, and the quantities Fn\qi, 
represent the standard PT interaction kernels, symmetrized in all of their arguments. We may 
obtain the time rate of change of the fluctuation by simply diff'erentiating the above, 

oo 

5(k, a) = f{a)H{a) ^ n[D{a)r6niK ao) , (4) 

n=l 

where f{a) = dlogD{a)/dloga. Inserting Eqs ([3]) and ([H) into Eq. ([2]), we find 

^k,a) = -n^o [^) ^ E [1 - ^/(«)] [D{arSn{k,ao) . (5) 

n=l 

Consider linear theory, n = 1: for the case O = 1, /(a) = 1 and potentials do not change with 
time and there is no ISW effect. However, for LCDM /(a) ~ and so 1 — /(a) > 0. Thus 
<^(k, a) > for an overdensity, and potentials decay at late times in linear theory. Consider 
now the nonlinear theory, n > 1: we immediately see that there are critical times in the LCDM 
model and that these are dictated by the sign of the bracket [1 — n/(a)]. For n = 2 we have 
1 — 2/(a) < for a < ajis, and 1 — 2/ (a) > for a > orsj and we call ors the Rees-Sciama 
time. For a < a^s nonlinear evolution enhances the growth of potential wells, but for a > crs 
it augments their decay. Thus it is theoretically possible to boost the late-time ISW effect. 



3 The ISW effect in A^-body simulations 

We use the Ziirich Horizon, "zHDRIZON", simulations to study the ISW. These are a large ensem- 
ble of pure CDM A^-body simulations (A^sim = 30). In this study we use the first 8 simulations, 
since these have 11 snapshots logarithmically spaced in the expansion factor from z = 1 to 
0, giving good time sampling. Each numerical simulation was performed using the publicly 
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Figure 1: Evolution of in a slab of thickness Ax = 100 h ^Mpc. Left z = 10 and right z = 0. 



available Gadget-2 code 10], and followed the nonlinear evolution under gravity of = 750^ 
equal mass particles in a cube of length L = 1500/i~^Mpc. The cosmological parameters are: 
{QmO = 0.25, = 0.75, ag = 0.8, Ug = 1.0, h = 0.72}, where these are: density parameters 
in matter and vacuum energy; power spectrum normalization and spectral index; Hubble pa- 
rameter. Dark matter halo catalogs were generated using the FoF algorithm, and the minimum 
number of particles for which an object was considered bound was 30. This gave haloes with 
M > 1.5 X lO^^MQ/h [for more details see I 111]. 

Considering again Eq. ([2]) we see that in order to measure the ISW in simulations we need 
to be able to estimate 6 and 5. The first is straightforward. The second can be obtained from 
the continuity equation: j(x) = —V • [1 + 5(x)] v(x)/a = —V • p(x)/o. In Fourier space this 
becomes, 6k = i^- Pfc/«, which may easily be computed [seelH, for details]. Figure [D presents a 
visual representation of the ISW effect in the simulations. 



4 Results: CMB-LSS angular cross-povi^er spectrum 



Under the Limber approximation the CMB-LSS angular cross-power spectrum [C^^'^) is [l|: 
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where P^^ = (5h(k)5^(— k)) is the 3D cross-power spectrum between halo density fluctuations 
5h and In Eq. ([6]), we have included the weight function 11^^, which for a mass-selected survey 
of clusters would have the form: 
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where is the comoving angular diameter distance and n(M, x) the cluster mass function. 
^TOTiXiiXj) = I^i dx^'^D\{x) Su^i^'^^'^i^'^x) 1 is the total number of clusters above mass 
-^min- The redshift shell is selected using the top-hat function: = [©(x ~ Xi) ~ ©(x ~ Xj)]-, 

with the Heaviside step function. 

Figure [2] presents the results for Cf^ between the CMB and group scale dark matter haloes 
(left panel) and cluster mass haloes (right panel). In each case we show the results for 5 narrow 
bins in redshift, with Az = 0.2. We find that on scales / < 100 the departures from linear theory 
are small < 10%, and are characterized by a small amplification of the signal, followed by a 
strong suppression. The departures appear small when compared to the cosmic variance, which 
is dominated by the Cf^ spectrum. 

We next investigated the S/M for and found good agreement with linear theory ex- 
pectations: the presence of bias effectively cancels out in the S/M expression and leads to 
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Figure 2: Angular cross-power spectrum of ISW effect and haloes as a function of spherical harmonic 
multipole L Left panel: results for group-scale dark matter haloes. Right panel: results for the cluster-scale 
haloes. In each panel we show results for 5 equally spaced bins in redshift over the range: z = [0.0, 1.0]. 
The predictions are differentiated by line thickness: thick lines - low redshift; thin lines - high. The 
line styles denote: linear theory - solid green line; nonlinear PT - red dash line; bi-cubic spline fit to 
the simulation data - blue triple-dot dash line. Top sections of each panel give the absolute power, and 
the lower sections the ratio with respect to linear theory. The shaded regions represents the l-cr error 
domains per multipole of the linear cross-spectra, where the central redshifts z E {0.1,0.3,0.5,0.7,0.9}, 
correspond to the colors (red, green, blue, cyan, magenta). 



negligible changes in the cross-correlation detectability. In fact, through the increased Poisson 
noise of the biased sample, there was a small reduction in the S/Af, relative to that for the dark 
matter-CMB cross-correlation. 

5 Conclusions 

We therefore conclude that the current power spectrum analyses of 0] and @] are not affected 
by nonlinear density evolution or scale-dependent bias and that these do not influence the 
detectability of the ISW-LSS cross-correlation. We are also led to believe that the results from 
0] are unlikely to be explained by nonlinear biasing effects in the LCDM model. Whether there 
remain systematic errors in the CMB data associated with point sources or whether this is an 
exciting detection of new physics remains an open question. 
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